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ABSTRACT 

Molecular-dynamics  simulations  of  activated  electron-transfer  (ET) 
reactions  in  methanol  have  been  undertaken  in  order  to  explore  the  physical 
nature  of  the  solvent  dynamics  coupled  to  ET  barrier  crossing,  and  to  probe  some 
underlying  reasons  for  the  facile  kinetics  observed  in  this  solvent.  The 
reactant  is  modeled  by  a  pair  of  Lennard-Jones  (LJ)  spheres  in  contact,  of 
varying  diameter  (4  or  5  A)  ,  and  containing  a  univalent  charge  (cation  or  anion) 
on  cne  site  so  to  probe  possible  effects  of  the  ionic  charge  sign.  Following 
equilibration,  the  collective  solvent  response  to  a  sudden  charge  transfer 
between  the  spherical  sites  is  followed,  and  described  in  terms  of  the  response 
function  C(t),  describing  the  difference  in  the  solvent-induced  electrostatic 
potential  between  the  initial  and  final  solute  states.  In  ail  cases,  the  C(t) 
curves  exhibit  a  very  rapid  (50-100  fs)  initial  decay  component  associated  with 
hydroxyl  inertial  motion,  followed  by  components  arising  from  hydrogen-bond 
librational  and  diffusive  motions.  Interestingly,  the  dynamics  and  relative 
importance  of  these  relaxation  modes  are  dependent  on  the  charge  sign  as  well  as 
size  of  the  solute  pair.  The  molecular-level  origins  of  these  sensitivities  are 
explored  by  examining  time-dependent  radial  distribution  functions,  which 
implicate  the  dominance  of  short-range  solvation  in  the  rapid  relaxation 
d3mamics.  In  particular,  the  initial  very  rapid  C(t)  component  for  the  smaller 


anion-neutral  reactant  pair  is  seen  to  arise  chiefly  from  dissipation  of  the 
hydroxyl  solvent  polarization  around  the  newly  formed  neutral  site  following  ET, 


being  accompanied  by  a  slower  build  up  of  solvent  structuring  around  the  adjacent 
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anionic  site.  The  simulated  ET  reorganization  energies  are  also  shown  to  be  TAB 
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dependent  upon  the  reactant  size  and  charge  type.  Some  more  general  implications 


of  these  and  other  MD  simulation  results  to  the  elucidation  of  dynamical  solvent 


effects  in  activated  ET  processes  are  also  noted. 
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Rates  of  chemical  reactions  involving  charge  redistribution  are  often 
influenced  strongly  by  polar  solvents.  Over  the  past  decade  there  has  been  a 
rapid  development  in  our  theoretical  and  experimental  understanding  of  the  role 
of  solvent  dynamics  on  such  reaction  rates.  These  studies  have  been  concerned 
with  time-dependent  fluorescence  Stokes  shifts  (TDFS),  ionic  and  dipolar 
solvation  dynamics [ 1-6 ] ,  and,  of  particular  interest  to  us  here,  activated 
electron-transfer  (ET)  reactions [ 7 , 8 ] .  Progress  in  these  closely  related  areas 
has  been  the  subject  of  a  number  of  informed  review  articles,  including  those 
cited  above.  The  dynamics  of  solvent  reorganization,  which  is  expected  to 
influence  adiabatic  ET  races,  ia  linked  intimately  to  the  response  of  the  solvent 
to  an  instantaneous  change  in  the  solute  charge  distribution  observed  in  TDFS 
experiments . 

While  often  beset  with  difficulties  in  separating  contributions  to  the 
observed  ET  rates  from  dynamics  and  activation  energetics,  there  is  nevertheless 
firm  experimental  evidence  that  overdamped  solvent  relaxation  can  dominate  the 
adiabatic  barrier-crossing  frequency,  i/„,  in  some  cases[7,8].  The  situation  is 
apparently  most  straightforward  in  so-called  Debye  solvents,  for  which  the 
uniform-phase  dielectric  relaxation  is  described  approximately  by  a  single 
exponential  decay.  For  reactions  in  such  media  where  the  activation  barrier  is 
also  associated  chiefly  with  solvent  reorganization,  is  found  to  correlate 
roughly  with  the  inverse  longitudinal  relaxation  time,  as  anticipated  on 
the  basis  of  the  simple  dielectric-continuum  picture(8].  A  notably  different 
situation  may  apply  at  least  to  reaction  dynamics  in  non-Debye  liquids, 
especially  alcohols,  that  feature  additional  shorter-time  components  in  the 
dielectric  loss  spectrvun.  Barrier-crossing  frequencies  that  are  markedly 
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accelerated  above  have  been  observed  for  a  number  of  outer-sphere  ET 
reactions  in  primary  alcohols ( 9-16 ] .  Part  of  the  often  large  (-10  fold)  rate 
accelerations  can  be  accounted  for [15)  In  terms  of  a  continuum-based  formulation 
due  to  Hynes [17].  A  subsequent  comparison  between  TDFS  data  obtained  for 
coumarins  and  ET  dynamics  extracted  from  kinetic  and  optical  data  for  metallocene 
self-exchanges  highlights  the  importance  of  rapid  (<  1  ps)  solvent  relaxation 
components  to  the  barrier-crossing  dynamics  in  methanol [ 7 , 16 ] . 

Such  documented  rate  effects  provide  additional  motivation  for  elucidating 
the  molecular-level  nature  of  the  remarkably  facile  dynamics  that  are  partly 
responsible  for  the  ET  reaction  rates  in  this  and  other  solvents.  To  this  end, 
it  is  necessary  to  proceed  beyond  dielectric-continuum  models  in  order  to  explore 
the  role  of  short-range  solvation  in  the  reaction  dynamics.  Theoretical 
approaches  to  this  problem  have  included  the  use  of  integral  equation  methods, 
such  as  the  mean  spherical  approximation  (MSA),  and  generalized  continuum  tactics 
which  attempt  to  account  for  the  wavevector  (k)-  as  well  as  frequency-dependence 
of  the  solvent  dielectric  permlttlvlty[3,4,18,19]  .  The  former,  however,  suffer 
from  an  inadequate  account  of  solvent-solvent  interactions,  while  the  latter  face 
the  difficulty  of  deducing  suitable  k-dependent  dielectric  parameters. 

Recently,  several  authors  have  explored  these  and  related  solvation  issues 
by  means  of  molecular  dynamics  (MD)  simulations [4 , 20-30] .  While  the  simulated 
systems  inevitably  provide  only  incomplete  models  of  the  actual  condensed-phase 
environment,  they  can  aid  considerably  our  appreciation  of  the  molecular- level 
details  of  nonequilibrium  as  well  as  equilibrium  solvation,  and  provide  a 
critical  test  of  analytic  theories  of  solvation  dynamics [4 , 5] .  A  significant 
common  finding  for  the  MD  studies  reported  to  date  is  that  fast,  non-diffusive , 
relaxation  mechanisms,  primarily  inertial  dynamics  and,  in  hydrogen-bonding 
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solvents,  librational  dynamics  as  well,  can  constitute  major  contributors  to  the 
solvent  response  to  charge  transfer.  Methanol  Is  the  only  solvent  studied  so  far 
for  which  these  fast  relaxation  mechanisms  do  not  dominate  the  solvation 
dynamics,  with  overdamped  (diffusive)  motion  being  an  important  additional 
contributor.  The  results  of  Fonseca  and  Ladanyi[20]  on  the  relaxation  of  the 
methanol  solvent  in  response  to  dipole  creation  in  a  diatomic  solute  differ  from 
other  MD  results  on  solvation  dynamics  in  another  important  respect:  a 
significant  departure  of  the  nonequilibrium  solvation  response  from  the 
relaxation  predicted  by  the  linear-response  (UR)  approximation  was  found.  Other 
MD  studies  found  good  agreement  with  LR  for  comparable  magnitude  of  the  change 
in  solute-solvent  interactions (4] .  This  breakdown  of  LR  in  methanol  is 
apparently  related  to  the  fact  that  several  relaxation  timescales  contribute 
importantly  to  the  solvation  dynamics ( 20] .  Since  this  characteristic  is  shared 
by  other  alcohols  as  well  as  other  non-Debye  solvents,  the  validity  of  LR,  which 
is  commonly  presumed  in  theoretical  models  of  solvation  dynamics,  has  been  called 
into  question. 

The  present  paper  reports  MD  simulations  of  intermolecular  ET  reactions  in 
methanol,  the  reactant  pair  being  modeled  by  a  pair  of  juxtaposed  spherical 
solutes,  one  charged  and  the  other  uncharged.  The  role  of  solute-solvent 
interactions  is  explored  by  altering  both  the  solute  size  and  univalent  charge 
sign,  thereby  mimicking  possible  behavioral  differences  between  the  cation- 
neutral  and  anion-neutral  reactant  pairs  which  commonly  constitute  the  self¬ 
exchange  process  utilized  to  probe  solvent-dynamical  effects [8].  The  solvent 
reorganization  d)rnamics  coupled  to  the  activated  reaction  are  simulated  by 
transferring  an  electron  between  the  thermally  equilibrated  solute  sites  so  to 
exchange  the  charge ,  and  following  the  ensuing  temporal  changes  in  solvent 
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Structure  and  energetics.  This  situation  corresponds  to  that  encountered 
experimentally  in  photolnduced  electron  transfer;  there  are  nonetheless 
intrinsically  close  relationships  between  the  dynamics  and  energetics  of 
photolnduced  and  thermally  activated  ET  processes.  A  similar  MD  simulation 
procedure  was  followed  earlier  by  Chandler  and  coworkers  in  their  examination  of 
pg3+/2+  exchange  in  water[26].  It  is  nevertheless  worth  noting  that  the  present 
MD-ET  simulations  differ  from  the  majority  of  previous  MD  studies  which  are 
concerned  with  the  formation  (or  annihilation)  of  a  dipole  from  a  neutral  solute 
(or  solute  pair) .  The  latter  case  refers  specifically  to  TDFS  and  related 
ultrafast  experimental  probes  of  solvation  dynamics.  While  inherently  related, 
the  effect  of  the  net  (unipolar)  reactant  charge  upon  the  solvation  dynamics  and 
energetics,  of  central  interest  in  intermolecular  (and  most  intramolecular)  ET 
systems,  does  not  appear  in  the  dipolar  systems. 

The  present  simulations  reveal  dynamics  which  uniformly  contain 
contributions  from  three  distinct  relaxation  mechanisms  -  inertial,  hydrogen-bond 
librational,  and  diffusive  motions.  The  relative  contributions  of  these 
mechanisms  to  the  overall  dynamics  vary  according  to  the  solute  size  and  charge 
sign  as  well  as  timescale.  The  observed  differences  in  the  early-time  dynamics 
are  shown  to  be  related  to  dissimilarities  in  the  observed  equilibrium-solvation 
structures;  the  variations  in  the  djmamlcs  of  different-size  anions  can  be 
rationalized  further  in  terms  of  time-dependent  radial  distribution  functions. 
We  also  test  the  linear-response  approximation  for  charge  transfer  in  methanol 
and  find  that  it  is  not  strictly  valid.  Finally,  we  comment  upon  the  more 
general  significance  of  the  present  findings  to  the  molecular-level  understanding 
of  solvent  dynamical  effects  in  electron  transfer. 
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MODELS  AND  METHODS 

The  model  system  used  in  our  MB  simulations  consisted  of  255  methanol 
molecules  and  one  two-phase  solute.  Intermolccular  interactions  were  represented 
as  site-site  Lennard-Jones  (U)  +  Coulomb  pair  potentials,  with  parameters  for 
the  methanol  and  the  solutes  as  given  in  Table  I.  The  Lorentz-Berthelot 
combining  rules  -  (Oj,  +  a^)/2  and  were  used  to  obtain  Lennard- 
Jones  (U)  diameters  and  energies  for  interactions  between  pairs  of  unlike 
sites  a  and  /9.  The  model  for  methanol  used  in  this  work  was  developed  by 
Haughney  et  al[31];  it  is  a  rigid  bond-length  and  bond-angle  treatment  in  which 
the  methyl  (Me)  group  is  represented  by  a  single  interaction  center,  and  the 
other  two  sites  are  the  oxygen  (0)  and  hydroxyl  hydrogen  (H)  atoms.  The  physical 
properties  around  room  temperature [ 31-35] ,  including  the  wavevector  dielectric 
permittivity  tensor  [32]  as  well  as  some  aspects  of  solvation  dynamics [20]  in  this 
model  solvent  have  previously  been  studied.  The  reactant  partners  are 
represented  as  a  pair  of  spheres  -  one  neutral,  the  other  charged  -  both  of  which 
have  the  same  mass  and  Lennard-Jones  parameters  e,  and  a,.  The  intersite 
distance  is  held  fixed  at  a,,  i.e.,  the  two  spheres  are  essentially  in  contact. 
Charge-transfer  dynamics  data  were  obtained  for  solutes  of  two  different 
diameters  (o,  -  4  A  and  5  A)  and  charges  (+e  and  -e)  .  (The  examination  of  larger 
solutes  in  the  present  study  was  precluded  by  the  markedly  longer  computational 
times  arising  from  the  greater  number  of  solvent  molecules  that  thereby  need  to 
be  included.)  In  the  remainder  of  this  paper,  these  solutes  will  be  referred  to 
by  their  a,  values  and  charge  signs  (e.g.  4  A  cation). 

In  the  simulation,  the  256  molecules  were  placed  in  a  cubic  box  with  sides 
of  25.9  A  length.  Periodic  boundary  conditions  and  a  spherical  cut-off  at  one- 
half  of  the  box  length  were  applied  to  the  Lennard-Jones  interactions.  A 
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modified  version  of  the  Ewald  sum,  which  excludes  contributions  from  solutes 
outside  the  central  box,  was  used  to  treat  long-range  electrostatic  interactions. 
The  equations  of  motion  were  integrated  using  the  Verlet  leapfrog  algorlthm[36] 
and  SHAKE[37]  was  utilized  to  maintain  the  bond  length  and  angle  constraints. 
The  production  runs  were  carried  out  in  the  NVE  ensemble  at  an  average 
temperature  of  298  K,  using  U  fs  time  steps. 

Two  types  of  MD  simulation  were  performed.  The  response  of  Che  system  to 
an  instantaneous  (q,0)  (O.q)  electron  transfer  between  the  charged  and  the 
neutral  solute  site  was  simulated  directly  by  nonequilibrium  MD.  In  addition, 
MD  trajectory  data  were  generated  to  calculate  the  equilibrium  solvent  structure 
around  the  solutes  and  to  obtain  the  linear-response  approximation  to  the  charge- 
transfer  perturbation. 

Each  starting  configuration  for  a  nonequilibrium  trajectory  was  generated 
from  an  initial  equilibrated  state  by  the  following  procedure.  The  charge  on  the 
solute  was  turned  off  so  that  the  surrounding  solvent  was  able  to  move  about  more 
freely  than  in  the  presence  of  a  charged  solute  and  the  system  was  allowed  to 
evolve  for  250  time  steps,  each  of  8  fs  duration.  The  charge  was  then  reinstated 
and  the  system  allowed  to  come  to  equilibrium  for  UOO  to  1400  additional  8  fs 
time  steps.  Finally,  since  the  actual  nonequilibrium  trajectories  were  obtained 
using  a  4  fs  step  size,  each  of  the  starting  configurations  obtained  as  described 
above,  was  equilibrated  for  an  additional  20  steps,  each  of  4  fs  duration.  The 
temperature  of  the  system  remained  stable,  indicating  that  nothing  unexpected 
happened  as  a  result  of  the  change  of  time  step  size.  The  equilibration  times 
for  different  solute-solvent  systems  were  not  the  same;  the  5  A  cation  and  the 
5  A  anion  solutes  were  equilibrated  for  400  time  steps,  the  4  A  cation  for  600 
time  steps  and  the  4  A  anion  for  1400  time  steps.  Longer  equilibration  times  for 
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the  smaller  solutes  and  especially  for  the  small  anion  system  reflect  the  fact 
that  a  "tighter"  solvation  shell  structure  needs  to  be  constructed  after  the 
ionic  charge  is  generated,  with  accompanying  slower  dynamics  (vide  infra). 

The  response  of  the  system  to  the  (q,0)  -►  (0,q)  charge-transfer 
perturbation  was  observed  for  a  period  of  600  fs.  At  least  200  trajectories  were 
averaged  to  obtain  the  response  function  of  each  solute.  The  quantity  calculated 
in  these  nonequilibrium  simulations  is  the  response  function 

C(t)  -  [AE(t)  -  AE(«)]/[AE(0)  -  AET^]  (1) 

where  AE  -  E,  -  Eq,  Eg  and  Eq  are  the  electrostatic  potentials  due  to  the  solvent 
at  the  charged  (q)  and  neutral  (0)  solute  sites,  t  -  0  is  the  time  at  which  the 
charge-transfer  perturbation  is  applied,  t  -  ®  corresponds  to  the  equilibrated 
solute-solvent  system,  and  the  overbar  denotes  an  average  in  the  presence  of  the 
perturbation,  i.e.  over  the  nonequilibrium  trajectories. 

As  noted  below,  comparison  of  the  above  response  function,  C(t),  with  the 
expectation  of  LH  theory  invites  the  calculation  of  the  time-correlation  function 

A(t)  -  (5AE(t)  5AE(0)>/((5AE)2)  (2) 

where  5AE(t)  -  AE(t)-(AE)  is  the  fluctuation  in  AE  and  (...)  denotes  an 
equilibrium  ensemble  average.  Because  the  two  solute  sites  differ  only  in  their 
charge  and  in  view  of  the  symmetry  of  the  charge-transfer  perturbation,  AE(®)  - 
(AE)  and  AE(0)  -  -(AE) .  For  each  solute-solvent  system,  the  time  correlation, 
A(t) ,  was  calculated  from  at  least  5  equilibrium  runs  of  7500  4  fs  steps.  The 
starting  points  for  these  equilibrium  runs  were  chosen  from  the  starting 
configurations  for  the  nonequilibrium  trajectories,  which  were  continued  without 
applying  the  charge-transfer  perturbation. 
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RESULTS  AND  DISCDSSION 

The  primary  results  of  this  work,  the  time  (t)  evolution  of  the 
nonequilibrium  solute-solvent  response  functions,  C{t),  for  the  four  solute  cases 
considered  here,  4  and  5  A  diameter  spherical  pairs  containing  univalent  cations 
and  anions,  are  presented  in  Figure  1.  For  clarity,  the  anion-response  functions 
are  offset  by  0.2  upwards  on  the  y-axis  from  the  cation  functions.  All  the  C(t) 
curves  exhibit  hallmarks  of  three  relaxation  mechanisms:  (1)  Inertial  dynamics, 
giving  rise  to  a  Gaussian  decay  at  short  times;  (2)  hydrogen-bond  librations, 
commencing  at  intermediate  times,  after  the  initial  inertial  component;  (3) 
strongly  overdamped  (dif fusional)  relaxation  at  longer  times.  The  presence  and 
characteristics  of  these  relaxation  modes  were  sketched  by  Fonseca  and  Ladanyi 
in  their  MD  simulation  of  dipole  creation  in  methanol [20] .  These  are  also 
evident  in  dipolar  relaxation  of  pure  methanol  solvent  as  observed  in  recent  far- 
infrared  absorption  experiments [38]  and  in  HD  simulations [ 35 ] .  The  same 
relaxation  mechanisms  contribute  to  solvation  dynamics  in  water[22,26,27 j ,  but 
there  diffusional  relaxation  apparently  plays  a  much  less  Important  ''ole[20]. 

Interestingly,  the  relative  contributions  of  these  different  relaxation 
mechanisms  for  the  present  ET  systems  are  seen  to  differ  according  to  the  solute 
size  and  charge  sign  (Figure  1) .  Anticipating  a  similarity  in  the  longer-time 
relaxation  rate  for  all  solutes,  the  overall  solvent  relaxation  time  is  shorter 
for  cations,  more  noticeably  for  the  smaller  solutes,  than  for  anions  despite  the 
fact  that  the  initial  inertial-driven  decay  is  more  rapid  for  the  latter.  These 
differences  in  C(t)  for  ions  of  the  same  charge  and  different  size  are 
significant,  but  do  not  exhibit  a  simple  trend:  faster  relaxation  rates  are 
observed  for  smaller  cations  and  larger  anions.  In  broad  outline,  these  C(t) 
traces  are  reminiscent  of  the  MD  simulations  obtained  for  dipole  creation  in 
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methanol  by  Fonseca  and  Ladanyi(20].  As  might  be  expected,  the  latter  C(t) 
responses  appear  to  lie  somewhat  in  between  the  C(t)  curves  for  the  cation  and 
anion  systems.  The  present  data,  however,  show  the  presence  of  a  more  rapid 
overall  decay  to  lower  C(t)  values  (<  0.4-0. 5),  especially  for  the  smaller  cation 
reactant . 

For  ion-neutral  solute  pairs  of  the  same  size,  the  nature  of  simple  ion- 
dipole  solute-solvent  interactions  should  be  independent  of  the  sign  of  the  ionic 
charge.  Consequently,  the  observed  differences  in  the  present  simulations  signal 
the  importance  of  more  involved  solute-solvent  interactions,  manifesting 
themselves  in  the  structure  (i.e.  molecular  orientation,  etc.)  of  the  first  one 
or  two  solvation  shells  and  in  the  rearrangement  dynamics  within  these  shells  in 
response  to  the  electron  transfer.  In  order  to  gain  a  more  specific  molecular- 
level  understanding  of  the  nature  of  the  C(t)  relaxations  and  their  dependencies 
on  solute  size  and  charge  sign,  it  is  therefore  instructive  to  inspect  the 
temporal  evolution  of  the  radial  distribution  functions  following  electron 
transfer,  as  also  can  be  provided  by  the  MD  simulations. 

The  various  site-site  radial  distribution  functions,  g(r) ,  for  the  ri,  0, 
and  Me  sites  of  methanol  solvating  the  4  A  charged  solute  spheres  (cation  or 
anion,  as  indicated)  under  equilibrium  conditions  are  plotted  in  Figure  2A;  the 
corresponding  data  for  the  5  A  ions  are  shown  in  Figure  2B.  Figure  3A  provides 
similar  distribution  functions  for  the  4  A  cation  and  anion,  but  now  measured 
from  the  center  of  the  uncharged  reactant  partner  within  the  charged-neutral 
reactant  pair.  Figure  3B  contains  corresponding  information  for  the  5  A  solutes. 

The  most  striking  aspect  of  Figure  2  is  the  large  difference  in  the  degree 
of  solvent  structuring  around  the  cations  and  anions.  Because  the  hydroxyl 
hydrogen  has  no  LJ  potential  and  carries  a  substantial  charge  (0.43e), 
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equilibrium  structures  yielding  a  close  proximity  between  H  and  the  ani in  solute 
are  strongly  favored,  resulting  in  a  relatively  ordered  solvent  environment 
(Figure  2).  In  t!  e  case  of  cations,  the  attractive  coulombic  interaction 
involves  the  methanol  oxygen,  for  which  the  distance  of  closest  approach  is 
significantly  larger  due  to  U  repulsion,  so  that  the  solvent  structuring  is  more 
diffuse.  As  might  be  expected  from  the  C(t)  data,  the  location  of  the  nearest- 
neighbor  g(r)  peaks  are  also  dependent  upon  the  ion  size.  For  the  cations,  the 
stronger  coulombic  attraction  yields  a  higher  and  narrower  first  oxygen  g(r) 
peak,  g_o,  surrounding  the  4  A  ion,  but  the  solvent  structure  around  the  4  and 
5  A  ions  is  otherwise  quite  similar.  In  contrast,  the  anionic  solutes  display 
a  marked  ?  nsitivity  of  the  g(r)  solvent  structure  to  the  ion  size.  For  the  4 
A  ions,  the  nearest-neighbor  -H  peak  is  very  high  and  sharp,  indicating  the 
predomirance  of  oriented  methanols  with  the  hydroxyl  H  pointiiig  towards  the  ion 
and  the  o  and  Me  sites  away  from  it.  In  the  case  of  5  A  anions,  however,  other 
structures  are  also  quite  likely,  as  can  be  deduced  from  the  relatively  sharp 
second  g_g  and  g_o  peaks  and  from  the  fact  that  the  first  g-^,  peak  is  higher  than 
the  first  g_g  peak.  This  is  indicative  of  structures  in  which  the  methyl  group 
orients  preferentially  towards  the  ion,  with  the  hydroxyl  hydrogen  po.sitioned 
either  tevards  or  away  from  the  charge.  The  latter  structure  allows  hydrogen 
bonding  between  solvent  molecules  in  the  first  and  second  solvation  shells, 
leading  to  more  favorable  solvent-solvent  interactions. 

As  expected,  the  solvent  ordering  around  the  adjacent  neutral  site  is  less 
sensitive  to  the  solute  size  and  the  charge  sign  of  the  adjacent  ionic  site 
(Figure  3).  The  neutral  site  is  seen  to  interact  most  favorably  with  the  methyl 
group,  which  has  the  largest  LJ  energy  and  the  smallest  charge.  Thus  the  first- 
neighbor  Me  peak  is  uniformly  the  largest  and  in  the  case  of  the  5  A  solutes,  at 
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the  shortest  intersite  distance.  For  the  5  A  solutes,  the  structure  around  the 
neutral  is  largely  uninfluenced  by  the  charge  sign  of  the  juxtaposed  luuIc  solute 
(Figure  3B) .  Some  effects  of  the  presence  of  the  adjacent  ion,  however,  are  seen 
for  the  4  A  solutes  (Figure  3a)  .  They  are  weak  for  the  0  and  Me  centers,  but  the 
H  distribution  favors  considerably  shorter  intersite  separation  when  the  neutral 
is  adjacent  to  the  anion.  This  is  consistent  with  the  fact  that  a  short 
hydrogen-solute  distance  is  strongly  favored  for  the  neighboring  anion. 

These  variations  in  the  local  equilibrium  solvation  for  differing  solute 
pairs  can  help  rationalize  some  of  the  aforementioned  differences  in  the  early- 
time  dynamics  following  electron  transfer  (Figure  1).  In  all  cases,  the  inertial 
component  is  due  primarily  to  the  motion  of  free  hydroxyl  hydrogens [20] . 
Immediately  after  the  electron  transfer,  the  density  of  hydroxyl  hydrogens  is  the 
largest  in  the  vicinity  of  the  newly  formed  neutral  site  of  the  4  A  anionic 
solute.  These  hydrogens  are  free  to  rotate  around  the  methanol  C-0  bonds,  giving 
rise  to  the  larger  initial  decay  component  of  C(t)  for  the  smaller  anionic  solute 
(Figure  1).  Hydrogen-bond  libration,  which  corresponds  to  frustrated  (l.e.  non- 
free)  rotation  of  the  hydroxyl,  plays  a  more  important  role  for  cations  at  early 
times,  because  the  hydrogens  in  the  first  solvation  shell  are  repelled  by  the  ion 
and  are  more  likely  to  form  bonds  to  the  more  distant  solvent  molecules.  This 
coupling  to  the  solvent  hydrogen-bond  network  is  initially  stronger  for  5  A 
versus  4  A  cations,  accounting  for  the  earlier  and  larger-amplitude  librational 
oscillation  for  the  larger  ion  (Figure  1). 

While  the  equilibrium  structures  around  the  two  solute  sites  can  help 
understand  some  aspects  of  the  dynamics,  they  are  less  useful  at  intermediate 
times,  where  nonequilibrium  (time-dependent)  structural  information  is  more 
pertinent.  In  this  latter  time  regime,  a  significant  difference  is  observed  in 
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the  solvent  relaxation  rate  for  the  U  and  5  A  anions,  but  less  so  for  the 
corresponding  cations  (Figure  1)  .  We  therefore  concentrate  on  the  former  in  the 
impending  discussion.  Figure  1  shows  that  the  response  function  for  the  A  A 
anion  decays  rapidly  up  to  about  80  fs,  exhibits  a  slower  relaxation  rate  to  200 
f s ,  and  then  accelerates  [relative  to  the  corresponding  C(t)  for  .he  5  A  anion] 
during  the  next  150  fs.  It  is  of  particular  interest  to  unravel  the  molecular- 
level  factors  responsible  for  such  complex  C(t)  behavior. 

A  useful  pictorial  means  of  displaying  information  pertinent  to  this  issue 
is  in  the  form  of  radial  distribution  functions  surrounding  both  solute  sites 
gathered  as  a  function  of  time  following  electron  transfer.  Figures  AA  and  B 
show  such  data  for  the  hydroxyl  hydrogen  surrounding  the  A  A  anion-neutral  pair 
for  the  indicated  series  of  time  periods  after  electron  transfer;  A  and  B  refer 
to  the  functions  for  the  newly  created  anion  and  neutral  sites,  respectively. 
[Note  that  the  initial,  t  -  0,  g(r)  functions  are  omitted  from  the  bottom  of 
Figures  AA  and  B  since  these  are  identical  to  the  equilibrixan  (t  -  «)  functions 
also  given  in  these  figures;  thus  the  t  -  0  function  for  the  newly  formed  neutral 
site  is  identical  to  the  t  -  ®  function  for  the  anion  site  (Figure  AA) ,  and  vice 
versa.]  Except  for  the  t  -  ®  data  which  were  obtained  from  equilibrium 
simulations,  these  quantities  were  extracted  from  nonequilibrium  trajectory  data. 
Due  to  the  statistical  limitations,  each  nonequilibrium  trace  refers  to  an 
appropriately  short  (28  fs)  time  interval. 

Inspection  of  Figure  AA.  shows  that  the  sharp  pronounced  g_H(r)  maximum  at 
about  r  -  2.3  A  if-  developed  substantially  only  for  t  is  200  fs  around  the  newly 
formed  anion.  Comparison  of  the  t  -  «  trace  in  Figure  AA  with  the  early-time 
curves  in  Figure  AB  shows  that  the  hydroxyl  H  structuring  around  the  newly  formed 
neutral  site  is  attenuated  severely  even  for  t  -  80  fs.  the  g(r)  traces 
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exhibiting  only  a  diffuse  broad  structure  throughout  the  later  times.  The  form 
of  the  corresponding  hydroxyl  oxygen  g(r)  traces  largely  mirrors  this  behavior, 
although  the  changes  are  less  marked  at  early  times. 

One  can  therefore  identify,  at  least  qualitatively,  the  origins  of  the  C(t) 
morphology  for  the  4  A  anion  system.  The  Initial  facile  C(t)  decay  to  about  80 
fs  is  associated  primarily  with  rapid  destruction  of  the  polarized  solvent 
structure  around  the  newly  formed  neutral  site,  whereas  the  ensuing  slower 
relaxation  arises  chiefly  from  solvent  structure  developing  around  the  newly 
formed  anionic  center.  The  shorter  timescale  for  dissipating,  as  opposed  to 
creating,  the  polarized  solvation  sphere  may  be  rationalized  on  the  basis  of  the 
anticipated  case  of  dissociating  a  relatively  ordered  solvent  structure  as 
compared  to  forming  it,  i.e.  on  entropic  grounds. 

Such  time-dependent  g(r)  curves  can  also  account  in  broad  terms  for  the 
markedly  different  C(t)  morphologies  seen  for  the  5  A  versus  the  4  A  anion 
systems  (Figure  1)  .  Figures  5A  and  B  contain  analogous  data  as  in  Figures  4A  and 
B,  for  the  newly  anionic  and  neutral  sites,  respectively,  but  for  5  A  solutes. 
As  will  be  anticipated  from  the  prior  discussion,  the  temporal  changes  in  solvent 
polarization  are  distinctively  milder  for  the  5  A  as  compared  with  the  4  A 
solute.  Thus  although  the  short  100  fs)  and  longer-time  segments  of  the  C(t) 
trace  for  the  5  A  solute  can  again  be  identified  chiefly  with  structural  changes 
surrounding  the  newly  formed  neutral  and  anionic  solute,  respectively,  the 
temporal  variations  in  g(r)  are  distinctly  milder.  In  particular,  this  behavior 
accounts  in  simple  fashion  for  the  slower  C(t)  decay  for  the  5  A  versus  4  A  anion 
for  t  >  200  fs  (Figure  1) .  Thus  comparison  of  Figures  5A  and  4A  shows  that  the 
relative  development  of  the  g_H(r)  maximum  centered  at  ca  2.3  A  and  2.9  A  for  the 
4  A  and  5  A  anions,  respectively,  in  the  time  regime  200-600  fs  is  significantly 
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less  pronounced  for  the  larger  solute. 

Some  other  interesting,  although  more  subtle,  pieces  of  molecular-level 
dynamical  information  can  also  be  extracted  from  such  temporal  g(r}  curves.  For 
example,  Figure  4  shows  that  significant  polarization  of  the  second  solvation 
shell  around  the  newly  formed  anion  develops  early  -  it  is  already  reasonably 
well  defined  within  the  earliest  time  Interval  displayed  (64-92  fs)  -  but  that 
the  first  shell  is  still  largely  unstructured.  The  first  solvation  shell  is 
formed  subsequently,  partially  by  removing  molecules  from  the  second  shell,  as 
deduced  from  the  opposite  temporal  dependencies  of  the  appropriate  g(r)  maxima 
(Figure  4A) .  It  should  be  recognized,  of  course,  that  the  quantitative  link 
between  the  time-dependent  g(r)  traces  and  the  overall  response  function  is  a 
complex  one  since  the  latter  necessarily  involves  a  detailed  convolution  of 
intermolecular  forces.  Nevertheless,  such  g(r)  data  can  clearly  provide  some 
intriguing  insight  into  the  individual  molecular  motions  responsible  for  the 
solvent  dynamics . 

Theoretical  models  for  solvation  dynamics  usually  assume  that  the  linear 
response  (LR)  approximation  applies  to  this  process  [4]  .  Fonseca  and  Ladanyi[20]  , 
however,  have  shown  that  this  approximation  is  surprisingly  poor  for  dipole 
creation  in  methanol.  Here  we  test  its  validity  for  electron  transfer  in 
methanol.  Figure  6  compares  the  time-evolution  of  corresponding  C(t)  and  its  LR 
counterpart,  A(t) ,  for  4  A  ions.  Figure  7  is  the  counterpart  of  Figure  6  for  5 
A  ions.  Noticeable  and  even  marked  differences  are  seen  between  corresponding 
C(t)  and  A(t)  traces,  especially  for  the  anionic  systems  and  for  the  smaller  (4 
A)  solute.  These  disparities,  which  are  generally  diagnostic  of  a  breakdown  in 
the  LR  approximation[4,39] ,  would  appear  to  be  associated  with  the  specific 
solvation  factors  discussed  above[20].  Interestingly,  however,  another  common 
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hydrogen-bonded  solvent,  water,  yields  reasonable  accordance  to  the  LR 
approxiination[22 , 26-28] .  The  LR  approximation  also  holds  for  the  ET  energetics 
in  water,  a  finding  exemplified  clearly  in  the  near-perfect  quadratic  energy 
surface  for  the  system  as  deduced  from  MD  simulations [ 26 ]  .  The  likely 
reasons  for  the  marked  behavioral  difference  between  methanol  and  water  have  been 
discussed  in  ref.  20.  All  three  moments  of  inertia  of  water  are  small  and 
correspond  primarily  to  hydroxyl  hydrogen  motion.  Only  one  such  moment  of 
inertia  is  present  in  methanol.  Thus  fast  inertial  decay  and  0-H  bond  libration 
are  much  more  efficient  in  relaxing  the  energy  in  water  than  in  methanol. 

In  addition  to  solvation  dynamics,  the  nonequilibrivun  MD  simulations 
necessarily  also  provide  information  on  the  energetics  of  activated  electron 
transfer [40] ,  Specifically,  the  {AE(®))  values  mentioned  above  equal  the 
potential  difference  between  the  Initial  equilibrium  and  ET  excited  states,  so 
that  the  quantity  |q(AE(«))|  can  be  identified  with  the  so-called  reorganization 
energy.  A,  for  electron  transfer.  Table  II  summarizes  reorganization  energies 
determined  in  this  manner  for  the  present  systems  in  methanol,  A^.  Listed 
underneath  are  the  corresponding  quantities,  Ajx;,  estimated  from  the  conventional 
dielectric-continuum  model  due  to  Marcus  [41] .  (The  latter  was  obtained  from  the 
usual  two-sphere  model,  by  setting  the  internuclear  distance  equal  to  twice  the 
reactant  radius.)  Comparison  of  the  corresponding  Aj®  and  A^c  values  reveals  that 
they  are  in  reasonable  agreement  for  both  the  larger  (5  A)  solutes  and  the 
smaller  cation-neutral  pair,  yet  the  A^  value  for  the  4  A  anion-neutral  reactant 
is  markedly  larger  than  Apc.  Note  that  while  the  dielectric-continuum  model 
necessarily  predicts  identical  A  values  for  equal-sized  anionic  and  cationic 
redox  couples  in  a  given  solvent,  only  the  An,  values  for  the  larger  solutes  are 
comparable,  the  smaller  anion  solute  yielding  a  Af®  value  that  is  about  40% 


16 


larger  than  for  the  cationic  reactant  (Table  II). 

This  latter  sensitivity  of  A|®  to  the  charge  sign  of  the  ion-neutral 
reactant  pair  provides  another  manifestation  of  the  differences  in  the  short- 
range  solvation  of  the  anion  and  cation  reactants  discussed  above .  The  larger 
value  for  the  4  A  anion  is  unsurprising  given  the  marked  changes  in  short- 
range  solvation  attending  electron  transfer  for  this  system,  together  with  the 
breakdown  in  the  LR  approximation.  Nevertheless,  the  nearly  equal  values  for 
the  5  A  cation  and  anion  systems  is  more  unexpected  on  this  basis.  Given  that 
the  LR  approximation  is  noticeably  less  valid  for  the  5  A  anion  than  for  the 
cation  reactant  (Figure  7,  vide  supra),  the  free  energy-reaction  coordinates 
should  be  decidedly  nonparabolic,  so  that  disparate  A^.^}  values  might  be  expected 
for  these  systems.  It  is  worth  noting  that  the  effective  diameters  of  the 
majority  of  ET  reactants  commonly  utilized  to  probe  solvent  dynamical  effects  are 
somewhat  larger  than  5  A;  for  example,  a  «  7-8  A  for  simple  metallocenes [ 9-11 ] . 
Consequently,  the  reasonable  agreement  between  the  A^p  and  Ap^  values  observed  for 
the  5  A  reactants  (Table  I)  could  be  construed  as  lending  support  to  the  validity 
of  the  dielectric  continuum  model  in  this  solvent. 

Such  agreement  between  A^  and  Apc  is,  however,  probably  misleading.  While 
the  model  for  methanol  used  in  the  HD  simulations  accounts  for  the  complex 
intermolecular  forces  involved  in  the  solute-solvent  interactions ,  it  lacks  a 
description  of  the  electronic  polarizability,  instead  tacitly  setting  the  solvent 
optical  dielectric  constant,  top,  equal  to  unity.  It  is  well  known  that  an 
important  (or  even  predominant)  contributor  to  the  energetics  of  nonequilibrium 
solvent  polarization,  and  hence  to  A,  is  associated  with  the  solvent  electronic 
polarizability,  as  described  by  e^p  in  the  dielectric-continuum  limit[41,42] . 
Since  increasing  the  electronic  polarizability  stabilizes  the  ET  transition 
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state,  setting  c^p  -  1  should  yield  A  estimates  that  are  too  large.  (The 
magnitude  of  the  discrepancy  should  be  roughly  a  factor  of  two  in  polar  solvents, 
for  which  -  1.75  to  2.4[42].)  Inclusion  of  this  factor  should  diminish  A^d 
below  the  corresponding  Auc  estimates.  Indeed,  near-infrared  optical  measure¬ 
ments  of  A,  Eop,  for  mixed-valence  ferrocenium-ferrocene  systems  in  methanol 
yield  values  that  are  ca  20%  below  Aqc,  even  though  in  most  other  polar  solvents 
(with  the  notable  exception  of  water)  Eop  is  close  to  (within  5-10%  of)  Aocl43]. 
The  origins  of  these  disparities  are  unclear,  but  can  be  accounted  for 
semiquantitatively  by  an  "non-local"  electrostatic  treatment,  whereby  A  is 
diminished  in  such  hydrogen-bound  solvents  by  the  reduced  ability  of  the 
"structured"  polarization  to  respond  to  local  alterations  in  the  electric 
field[44] .  The  solvent  structuring  effects  observed  in  the  present  MD 
simulations  have  a  rather  different  physical  origin,  arising  from  short-range 
solute-solvent  interactions,  and  would  appear  to  increase  rather  than  diminish 
the  ET  barrier. 

Implications  For  ET  Solvent  Dynamics 

A  clear  feature  that  emerges  from  the  present  MD  simulations,  which  is  also 
seen  in  related  dipole-creation  and  charge-transfer  simulations  in  water {22,26- 
28]  and  acetonitrile [ 21]  ,  is  the  importance  of  a  notably  rapid  inertial  decay  to 
the  overall  relaxation  dynamics [4].  Indeed,  a  typical  finding  from  the  MD 
simulations  in  all  three  solvents  is  that  the  response  function  decays  to  below 
0.5  in  about  100  fs  or  less.  In  physical  terns,  then,  a  substantial  relaxation 
of  the  excess  energy  present  in  the  ET  transition  state,  l.e.  to  energies 
markedly  (;fc  several  kflT)  below  the  transition  state  in  the  products'  well,  can 
be  accomplished  extremely  rapidly  in  these  media.  This  rapid  relaxation  is  seen 
to  involve  the  first  one  or  two  solvation  shells  surrounding  the  solute.  The 
observation  of  especially  fast  dynamics  for  short-range  solvation  differs  from 
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the  well-known  Onsager  "snowball"  conjecture[45] ,  Slower  relaxation  dynamics  for 
"close-in"  solvent  are  also  predicted  by  the  nonequilibrium  MSA  treatment  based 
on  overdamped  relaxation[46] ,  and  have  been  observed  in  MD  simulations  under 
similar  conditions  to  the  MSA  model  solventI29].  Nevertheless,  rapid  short-time 
dynamics  can  be  predicted  under  some  conditions  from  analytic  theory,  as  well  as 
MD  simulations,  by  considering  partly  inertial  rather  than  purely  overdamped 
motions [47] . 

Given  that  there  is  clear  theoretical  evidence  that  higher-frequency 
relaxations  can  dominate  the  ET  barrier-crossing  frequency,  »/„,  when  they  make 
moderate  or  large  contributions  to  the  energetic  response [ 15 , 17 ,48 ] ,  one  is 
tempted  to  assert  that  the  reaction  dynamics  in  experimental  ET  systems  are  at 
least  partly  controlled  by  such  rapid  inertial  motions  rather  than  the  largely 
overdamped  dynamics  which  are  usually  considered  in  discussions  of  ET 
kinetics[8].  Of  course,  as  already  mentioned,  methanol  is  a  "non-Debye"  solvent, 
exhibiting  at  least  one  additional  higher-frequency  relaxation  in  the  dielectric- 
loss  spectrum[49]  .  This  is  apparently  manifested  as  a  rapid  solvation  relaxation 
component  (r,  -  1  ps)  in  ultrafast  TDFS  measurements  reported  by  Barbara  et  al 
in  methanol [16, 50 ] .  The  high-frequency  component  observed  in  the  present 
simulations  in  methanol,  r,  -  0.1  ps,  is  clearly  much  faster,  occurring  on  a 
timescale  below  that  readily  accessible  to  TDFS  experiments,  and  is  not 
associated  directly  with  the  "non-Debye"  properties  of  the  pure  methanol  solvent. 
Nevertheless,  the  longer-time  (i.e.  slower)  component  of  the  C(t)  traces  in 
Figure  1  exhibits  a  relaxation  time,  ca  0.5  ps,  which  is  compatible  with  the  fast 
component  observed  in  the  Barbara  et  al  measurements. 

*  It  is  worth  noting  that  the  relaxation  times  provided  Ey  the  present 
simulations  are  anticipated  to  be  about  twofold  too  short  as  a  result  of  the 
neglect  of  the  dielectric  function  of  the  solvent  electronic  polarizability  to 
the  overall  longitudinal  dynamics[4]. 
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A  key  unresolved  Issue  is  the  extent  to  which  the  remarkably  rapid 
solvation  dynamics  as  observed  by  MD  simulations  so  far  in  water,  acetonitrile, 
and  methanol,  are  present  in  other  media  known  (or  at  least  anticipated)  to  offer 
higher  solvent  friction.  The  notably  facile  barrier-crossing  frequencies  for 
adiabatic  electron-exchange  pro:esses  involving  suitable  redox  couples  such  as 
metallocenes  in  water,  acetonitrile,  and  methanol [10, 11, 15, 16]  are  certainly 
compatible  with  the  involvement  of  such  rapid  MD  relaxations ,  although  the  onset 
of  nonadiabiticity  anticipated  towards  higher  frequencies  may  limit  their 
Importance  except  in  ET  processes  featuring  unusually  strong  donor-acceptor 
electronic  coupling[16] .  Slower  ET  reaction  dynamics  are  observed  in  higher- 
friction  solvents,  such  as  benzonitrile  and  dimethylsulfide ,  so  to  yield  solvent- 
dependent  barrier-crossing  frequencies  varying  in  rough  accordance  with 
i.e.  as  expected  from  the  overdamped  dielectric-continuum  treatment [ 7,8 j .  This 
finding  by  itself  suggests,  and  has  been  conanonly  taken  to  indicate,  that 
overdamped  solvent  relaxation  plays  a  key  role  in  controlling  solvent  dynamical 
effects  on  activated  ET  processes.  However,  in  the  light  of  the  MD  simulations, 
an  alternative  (albeit  perhaps  less  palatable)  interpretation  is  that  the 
solvent-dependent  ET  dynamics  are  not  described  by  continuum-like  overdamped 
motion  at  all,  but  rather  by  the  faster  dynamics  of  short-range  solvation  which 
may  happen  to  correlate  roughly  with  Some  related  issues  along  these  lines 

will  be  pursued  elsewhere [51]  .  The  exploration  of  this  possibility  would  be  one 
motivation  for  pursuing  MD  simulations,  and  also  appropriate  analytic  theoretical 
treatments,  in  ostensibly  higher-friction  solvent  media. 
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TABLE  I. 

Potential  Parameters  for  the  Solvent  and  Solute 


Solvent® 


Site 

q/e^ 

(eAB)/K‘' 

ag/Ad 

mass/amu® 

Me 

0.297 

91.2 

3.861 

15.024 

0 

-0.728 

87.9 

3.081 

16.000 

H 

0.431 

0 

0 

1.008 

Solute^ 

Type 

q/e* 

(es/A)/KC 

o/k^ 

mass/amu® 

4  A  cation 

+1,0 

150 

4.0 

60.0 

4  A  anion 

-1.0 

150 

5.0 

60.0 

5  A  cation 

+1.0 

150 

4.0 

60.0 

5  A  anion 

-1.0 

150 

5.0 

60.0 

a)  geometry: 

^•-0  “ 

1.43  A;  Ro-h  -  0.945  A;  Angle„,_o-H  - 

108'’53'  . 

b)  Net  charge  on  site  indicated  in  left-hand  column. 

c)  Lennard-Jones  energy  parameter. 

d)  Lennard-Jones  diameter. 

e)  Site  mass. 

f)  Solute  is  a  pair  of  spheres  at  a  center-to-center  distance  of  05,  i.e. 
approximately  in  contact,  one  containing  the  charge  as  noted  in  left-hand 
column. 


Ih 


TABLE  II. 

Calculated  Electron-Transfer  Reorganization  Energies  (kcal  moI~^) 


Solute^ 

A  A  anion 

4  A  cation 

5  A  anion 

5  A  cation 

58.6 

42.5 

30.2 

31.4 

AA.5 

44.5 

35.5 

35.5 

a)  Solute  is  charged-neutral  spherical  pair  in  contact,  with  diameters  as 
noted . 

h)  Reorganization  energy  as  estimated  by  MD  simulation  (see  text). 


c)  Reorganization  energy  as  estimated  from  Marcus  dielectric-continuum  theory, 
using  solvent  parameters  as  in  ref.  9. 
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FIGURE  CAPTIONS 
Figure  I 

Response  function,  C(t),  traces  versus  time  following  instantaneous  symmetrical 
electron  transfer  within  ion-neutral  solute  pairs  in  methanol.  The  full  and 
dashed  lines  refer  to  4  and  5  A  diameter  solutes,  respectively.  Note  that  the 
traces  for  the  anion-neutral  pairs  are  offset  by  0.2  upwards  on  the  y-axis  for 
clarity. 

Figure  2 

A)  Equilibrium  radial  distribution  functions  for  the  H  (full  line) ,  0  (dashed 
line),  and  Me  (dotted  line)  sites,  as  indicated,  surrounding  the  charged 
sites  in  the  4  A  cation-neutral  and  anion-neutral  sites  (upper  and  lower 
panels,  respectively). 

B)  As  for  A),  but  for  5  A  solutes. 

Figure  3 

A)  As  for  Figure  2A,  but  for  neutral  4  A  site. 

B)  As  for  Figure  2B,  but  for  neutral  5  A  site. 

Figure  4 

A)  Time-dependent  radial  distribution  functions  for  hydroxyl  H  around  newly 
created  anion  site  following  electron  transfer  in  4  A  anion-neutral  solute. 
Traces  corresponding  to  increasing  times  are  offset  progressively  by  2.0 
upwards  along  the  y-axis . 

B)  As  for  A),  but  for  newly  formed  neutral  site. 

Figure  5 

A)  As  in  Figure  4A,  but  for  5  A  anion  site. 

B)  As  in  Figure  4B,  but  for  5  A  neutral  site. 
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Figure  6 

Comparison  of  the  non-equilibrium  response  function  C(t)  and  the  corresponding 
time-correlation  function  A(t)  vs  time  for  4  A  ion-neutral  pairs  in  methanol. 
Anion-neutral  curves  (full  lines);  cation-neutral  curves  (dashed  lines). 


Figure  7 

As  for  Figure  6 ,  but  for  5  A  ion-neutral  solute  pairs . 
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